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I outline the basic strategies for the computation of charm and bottom quark masses by means of lattice QCD, 
where particular emphasis is placed on the non-perturbative renormalization of the effective theory for the b-quark 
in heavy-light systems. A few selected results in the quenched approximation are reviewed, and the current status 
of extending these calculations to QCD with dynamical quarks is summarized. 



1. Introduction 

The fundamental parameters of QCD are the 
strong coupling constant and the quark masses. 
In particular, high-precision tests of the SM crit- 
ically depend on the masses of the charm and 
bottom quarks, and the non-perturbative nature 
of lattice QCD, which allows to directly compute 
hadronic observables to match experimental in- 
puts, is especially suited to determine them. 

Their ab initio lattice calculation faces system- 
atic errors, such as cutoff and finite-volume ef- 
fects, too large dynamical light quark masses and 
the omission of some (or, in quenched QCD, even 
all) sea quark flavours. Among these, the prob- 
lem of discretization errors proportional to pow- 
ers of the bare quark mass is of utmost relevance 
in case of systems containing a heavy quark, when 
a Wilson-like fermion action is employed. Requir- 
ing am qj h <C 1, it is then obvious that simulat- 
ing a b-quark at its physical mass of w 5 GeV in 
a space-time volume of i 4 = (2fm) 4 would de- 
mand 3> 50 lattice points per direction, an effort 
beyond the computing resources available today. 

In this situation, an attractive possibility is to 
recourse to an effective theory for the b-quark. 
Extracting physical predictions from it, however, 
involves a matching to QCD, which becomes a 
severe source of uncertainty when performed per- 
turbatively, for in the continuum limit results are 
power divergent. Before coming to the effective 
theory framework to deal with heavy-light sys- 



* Invited talk at the Workshop on e+e Collisions from $ 
to * (PHIPSI08), April 7-10, 2008, INFN Frascati, Italy. 



terns on the lattice (Section [3j) , let us consider the 
c-quark, where direct simulations are just doable. 

2. The charm quark's mass 

Physics in the region around the charm quark, 
which is by a factor of about 4 lighter than the 
b-quark, can be studied with conventional lattice 
QCD methods, provided that — assuming Wilson 
fcrmions here — the leading O(a) cutoff effects 
have been eliminated non-perturbatively and a 
range of lattice resolutions a -1 w (2 — 4) GeV 
is covered in the simulations. This is nicely il- 
lustrated by the benchmark calculation [1] of the 
mass of the c-quark in the quenched approxima- 
tion, which accounts for all systematic errors ex- 
cept for the neglection of dynamical quark. Fol- 
lowing the strategy and techniques used to com- 
pute m s [2|, and 0(am qiC ) effects oc (6a — bp) 
being removed [3], the renormalization group in- 
variant (RGI) c-quark mass was calculated on 4 
lattices with spacings a = (0.1 — 0.05) fm. In addi- 
tion to the definition of the (non-perturbatively) 
renormalizcd quark mass via the heavy-light 
PCAC relation, two further definitions were con- 
sidered, differing by 0(a 2 ) discretization errors. 
The agreement of all three definitions in the con- 
tinuum limit found in [T] thus provides a clear 
demonstration that on lattices of the chosen size 
(from 16 3 x 32 to 32 3 x 64) one can control the dis- 
cretization errors carefully by an extrapolation. 

This is confirmed by including a further gauge 
field ensemble in a volume of 48 3 x 96 with an 
even finer resolution of a = 0.03 fm, which was 
first generated in the context of [3] to investigate, 
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Figure 1. Continuum extrapolation of the RGI 
c-mass for Nf = 0, covering a = (0.1 — 0.03) frn; 
see [5j for details. The length scale r = 0.5 fm is 
derived from the force between static quarks. 

within quenched QCD, the size of cutoff effects 
in the charm sector thoroughly also for other ob- 
servables such as the D s -meson decay constant. 
Figure [1] shows an update of the continuum ex- 
trapolation of the RGI charm mass after a signif- 
icant increase of the statistics compared to |1I4) . 
particularly for the smallest a 5j. The final re- 
sult of this analysis, m c MS (m c ) = 1268(24) MeV, 
is compatible with the earlier result of ref. [1] . 

However, from the figure it is also evident that 
a controlled assessment of the cutoff effects will 
not be possible for much heavier quarks. Hence, 
for b-quarks other approaches must be devised, 
and one will be discussed in the next section. 

Charm-quark mass computations in QCD with 
dynamical flavours are not yet in the stage such 
that a solid continuum limit extrapolation can 
be performed [617] . Since non-perturbative es- 
timates for all the necessary improvement coef- 
ficients and renormalization factors in the two- 
flavour Wilson theory have become available now 
[8)9)10)11] . Nf = 2 calculations of masses and ma- 
trix elements in the charm sector along the lines 
of the quenched studies are being started [12] . 

For other recent determinations of m c using 
QCD sum rules and current-current correlators in 
lattice and continuum QCD, see refs. [13114] . 

3. Non-perturbative HQET 

Heavy Quark Effective Theory (HQET) at zero 
velocity on the lattice [15] offers a reliable solu- 
tion to the problem of dealing with the two dis- 
parate intrinsic scales encountered in heavy-light 
systems involving the b-quark, i.e., the lattice 
spacing a, which has to be much smaller than 



1 /nib to allow for a fine enough resolution of the 
states in question, and the linear extent L of the 
lattice volume, which has to be large enough for 
finite-size effects to be under control. 

Since the heavy quark mass (m\,) is much larger 
than the other scales such as its 3-momentum or 
Aqcd ~ 500 MeV, HQET relies upon a system- 
atic expansion of the QCD action and correlation 
functions in inverse powers of the heavy quark 
mass around the static limit (rrib — > oc). The 
lattice HQET action, Shqet, at 0(l/mb) reads: 

^Y/x^h {A) + 5m - cj kin D 2 - cj spin erB} ip h , 

with Vh satisfying P+ip h = vfa, P + = 
and the parameters Wkin and U! sp - m being formally 
0(l/mb). At leading order (static limit), where 
the heavy quark acts only as a static colour source 
and the light quarks are independent of the heavy 
quark's flavour and spin, the theory is expected to 
have ~ 10% precision, while this reduces to ~ 1% 
at 0(l/mb) representing the interactions due to 
motion and spin of the heavy quark. As crucial 
advantage (e.g., over NRQCD), HQET treats the 
1 /mb-corrections to the static theory as space- 
time insertions in correlations functions, viz. 

(O) = (O)stat + 

a 4 ^ {^kin (OOkin (X) ) st at + W spin (OO sp i n (x)) stat } 
x 

for multi-local fields O, where (0) s tat denotes 
the expectation value in the static approximation 
and Okin and O sp i n are given by 7/> h D 2 i/>h and 
i/> h <TB-0h- In this way, HQET at a given order is 
(power-counting) renormalizable and its contin- 
uum limit well defined, once the mass countert- 
erm 8m and the coefficients Wkin and w sp in are 
fixed non-perturbatively by a matching to QCD. 

Still, for lattice HQET and its numerical appli- 
cations to lead to precise results with controlled 
systematic errors in practice, two shortcomings 
had to be left behind first. 

1. The exponential growth of the noise-to- 
signal ratio in static-light correlation func- 
tions, which can be overcome by a clever 
modification of the Eichten-Hill discretiza- 
tion of the static action [TJ] . 
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2. As in HQET mixings among operators of 
different dimensions occur, the power- 
divergent additive mass renormalization 
dm r~j g^/a already affects its leading or- 
der. Unless HQET is renormalized non- 
perturbatively jf7j . this divergence — and 
further ones ~ g$ /a? arising at 0(l/mb) 
— imply that the continuum limit does not 
exist owing to a remainder, which, at any fi- 
nite perturbative order [18119] . diverges as 
a — > 0. A general solution to this theoret- 
ically serious problem was worked out and 
implemented for a determination of the b- 
quark's mass in the static and quenched ap- 
proximations as a test case in ref. [2D]. It 
is based on a non-perturbative matching of 
HQET and QCD in finite volume. 



In Section 13.11 1 review the HQET computa- 
tion of the mass of the b-quark including the 
0(l/m b ) terms for N t = 21J, while Section[3J2] 
briefly summarizes the present status of the on- 
going project to extend this to the more realistic 
case of two- flavour QCD [22] . A first application 
of the strategy of [2DJ to the computation of the 
B-meson decay constant can be found in [23] . 

A promising, not unrelated approach to lattice 
heavy-light systems derives from the step scaling 
method in the relativistic theory proposed in [24] , 
where extrapolations in the heavy quark mass of 
step scaling functions (SSFs) as finite-size effects 
of proper observables, combined with SSFs calcu- 
lated in HQET, turn into safer interpolations [25] . 

3.1. The b-quark mass in HQET at Of — ) 

As pointed out in 26], let us first note that 
in order not to spoil the asymptotic convergence 
of the series, the matching must be done non- 
pcrturbatively — at least for the leading, static 
piece — as soon as the l/m b -corrections are in- 
cluded, since as — ► oo the perturbative trun- 
cation error from the matching coefficient of the 
static term becomes much larger than the power 
corrections ~ Aqcd /wib of the HQET expansion. 

In the framework introduced in ref. [20] , match- 
ing and renormalization are performed simultane- 
ously and non-perturbatively. The general strat- 
egy, illustrated in figure [2J can be explained as 
follows. Starting from a finite volume with £j » 
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Figure 2. Idea of lattice HQET computations 
via a non-perturbative determination of HQET 
parameters from small-volume QCD simulations. 
Arrows indicate steps to be repeated at smaller a 
to reach a continuum limit. (Drawing from [261.) 

0.5 fm, one chooses lattice spacings a sufficiently 
smaller than 1/mb such that the b-quark prop- 
agates correctly up to controllable discretization 
errors of order a? . The relation between the RGI 
and the bare mass in QCD being known, suitable 
finite- volume observables $fc(Li,Mb) can be cal- 
culated as a function of the RGI b-quark mass, 
Mb, and extrapolated to the continuum limit. 
Next, the power-divergent subtractions are per- 
formed non-perturbatively by a set of match- 
ing conditions, in which the results obtained for 
$fc are equated to their representation in HQET 
(r.h.s. of figure [2]). At the same physical value 
of L\ but for resolutions L\/a = O(10), the 
previously computed heavy-quark mass depen- 
dence of $fc(ii,Mb) in finite- volume QCD may 
be exploited to determine the bare parameters of 
HQET for a » (0.025 - 0.05) fm. In order to 
evolve the HQET observables to large volumes, 
where contact with experiments can be made, one 
also computes them at these lattice spacings in a 
larger volume, L2 = 2L\. The resulting relation 
between <&/-(Lx) and $fe(i 2 ) is encoded in asso- 
ciated SSFs <7fc, as indicated in figure [2] Finally, 
the knowledge of $k{L2, Mb) and employing res- 
olutions Li I a = O(10) fixes the bare parameters 
of the effective theory for a » (0.05 — 0.1) fm 
so that a connection to lattice spacings is estab- 
lished, where large- volume observables, such as 
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the B-meson mass or decay constant, can be cal- 
culated (l.h.s. of figure [5]). This sequence of steps 
yields an expression of me (the physical input) as 
a function of My, via the quark mass dependence 
of $k(Li, Mb), which eventually is inverted to ar- 
rive at the desired value of the RGI b-mass within 
HQET. The whole construction is such that the 
continuum limit can be taken for all pieces. 

More specifically, upon restricting to spin- 
averaged quantities to get rid of the contribu- 
tions proportional to Lu sp - m , the task is to fix 
5m and Wkin non-perturbatively by performing 
a matching to QCD. For sensible definitions of 
the required matching observables, <&i and $2, 
we work with the Schrodinger functional (SF), 
i.e. QCD with Dirichlet boundary conditions in 
time and periodic ones in space (up to a phase 9 
for the fermions): $^ CD (I,mh) exploits the sen- 
sitivity of SF correlators to 9 and $2 D (L, mi,) = 
LTi(L, TOh), where Ti is a finite- volume effective 
energy. When expanded in HQElH, $f QET (L) 
is given by Wkin times a quantity defined in the 
effective theory (called RY n (L,9,9')), whereas 
$ 2 1C ^ ET (L) is a function of oj^ n and rribarc = 
dm + mi, involving two other HQET quantities, 
Tf at (i) and lf n (L). According to the strat- 
egy sketched above, by equating $j? CD (Li, mh) 
and <E>£ ET (Li) one can determine the bare pa- 
rameters mbaro and Wkin as functions of m-h at 
the lattice spacings belonging to the volume L\ . 
To use the spin- aver aged B-meson mass, m", as 
phcnomenological input, the <5>k are evolved to 
larger volumes through proper SSFs [21j . where 
the resulting $^ ET (2Li, mh) still carry the de- 
pendence on mh inherited from the matching to 
QCD in L\. After 2 evolution steps (and taking 
continuum limits), linear extents of > 1.5 fm are 
reached, and mbaro and Wkin, expressed in terms 
of SSFs, $^ CD (Li,m h ) as well as ^kin^ pstat and 

r^ m , are obtained — again as functions of m;,. 
Now, the b-quark mass is extracted by solving 

m av = £Stat +£Jk . n ( mh )£*in + (1) 

for m h , with E stat = limr^^ If at and E kin = 
-(B|a 3 J2 Z Okin(0, 2)|B)stat- All quantities enter- 



ing eq. (fTJ) having a continuum limit either in 
QCD or HQET implies that all power divergences 
have been subtracted non-perturbatively. 
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Figure 3. Graphical solution of eq. |3p in 
the quenched approximation [211. The quan- 
tity used in the hnite-volume matching step is 
$J CD (Li,M) = LxIx(L\,M), and z = L\M. 

In case of the leading-order, static approxima- 
tion, where only mbarc needs to be determined, 
the small- and large-volume matching conditions 



simplify to ri(Li,m h ) = r? at (Li) 



m ba 



and 



mff = £; stat + m ba re, respectively. To be able 
to replace mbare in the latter by the former, we 
bridge the volume gap in two steps by inserting 
a SSF cr m (ii) = 2Li[rf at (2Li) - rf at (Li)] and 
arrive at the master equation 



L\ [m 



\\ ^^stat pstat^j 



2 



Lili, (2) 



2 Here, 5m = in the action; its contribution is accounted 
for in the overall energy shift »rib arc in HQET versus QCD. 



where Ti ~ ri(Li, mj,) stems from Lf-QCD and 
any reference to bare parameters has finally dis- 
appeared. Its graphical solution is reproduced 
in figure [3] and yields M{f at = 6.806(79) GeV. 

For the details on the (technically more in- 
volved) inclusion of the sub-leading 1/mb-effects, 
which exploits the freedom in choices for the angle 
9 as well as an alternative set of matching observ- 
ables, I refer to ref. [21] . Here, I just quote their 

final result m^mb) = 4.347(48) GeV with the 
remark that, upon including the 1/mb-terms, dif- 
ferences among the static results w.r.t. the match- 
ing condition chosen are gone, which signals prac- 
tically negligible higher-order corrections. 

3.2. Status in two-flavour QCD 

The renormalization of HQET through the 
non-perturbative matching to Nf = 2 QCD in fi- 
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nite volume, to perform the power-divergent sub- 
tractions, is under way |22j . As an important pre- 
requisite, the non-perturbative relation between 
the RGI and subtracted bare heavy quark mass 
calculated in [101 1 lj enables to fix RGI heavy 
quark masses in the matching volume L\ : 

L\M = Z M {go)Z(g ) (1 + b m (g )am qth ) £im qih . 

The extent L\ is defined via a constant SF cou- 
pling, g 2 (L 1 /2) = 2.989, and the PCAC masses 
of the dynamical light quarks are tuned to zero. 

Figure U shows an example for the mass depen- 
dence of a finite-volume QCD observable in the 
continuum limit, which enters the matching step. 

L,r aB = L 1 (r ps +3r v )/4, 0=0.5 L,r av (z,e) 

13 |T= r n I 1 13 K I I I I I 1 

! E : I ': 

12 - • = : 

j — z-i5i 11 r 1 

10=, — « . ' z.13-: : : 

; ' " — • ; ,,-r \ 



0.001 0.002 (a/L,)* 6 9 10 11 13 13 14 15 L,M 

Figure 4. Preliminary continuum limit (left) 
and z-dependence (right) of the spin-averaged B- 
meson energy in hnite-volume QCD for N{ = 2. 

4. Outlook 

Thanks to technical and conceptual advances, 
most notably regarding non-perturbative O(a) 
improvement and renormalization in QCD and 
HQET, heavy-light physics has recently seen a 
significant progress. Since most current simula- 
tions control other systematics such as dynam- 
ical quark effects and the continuum extrapola- 
tion, more precise results for the c- and b-quark 
masses can be expected within the next years. 
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